Response Curves and Preimage Sequences 
of Two-Dimensional Cellular Automata 



Henryk Fuks and Andrew Skelton 

Department of Mathematics, Brock University, St. Catharines, Ontario, Canada. 



Abstract — We consider the problem of finding response 
curves for a class of binary two-dimensional cellular au- 
tomata with L-shaped neighbourhood. We show that the 
dependence of the density of ones after an arbitrary number 
of iterations, on the initial density of ones, can be calculated 
for a fairly large number of rules by considering preimage 
sets. We provide several examples and a summary of all 
known results. We consider a special case of initial density 
equal to 0.5 for other rules and compute explicitly the 
density of ones after n iterations of the rule. This analysis 
includes surjective rules, which in the case of L-shaped 
neighbourhood are all found to be permutive. We conclude 
with the observation that all rules for which preimage curves 
can be computed explicitly are either finite or asymptotic 
emulators of identity or shift. 
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1. Introduction 

Cellular automata (CA) can be viewed as computing 
devices, which take as an input some initial configuration. 
The CA rule is iterated a number of times, resulting in a final 
output configuration. In many practical problems, e.g., in 
mathematical modelling, one wants to know how a CA rule 
iterated over an initial configuration affects certain aggregate 
properties of the configuration, such as, for example, the 
density of ones. If we take a randomly generated initial 
configuration with a given density of ones, and iterate a given 
rule n times over this configuration, what is the density of 
ones in the resulting configuration? We want to know the 
"response curve", the density of the output as a function of 
the density of the input. 

Response curves appear in computational problems, and 
a classical example of such a problem in CA theory is the 
so-called density classification problem (DCP). If we denote 
the density of ones in the configuration at time n by P n (l), 
the DCP asks us to find a rule for which Poo(l) = 1 if 
P (l) > 1/2 and P 0O (l) = if P (l) < 1/2. Since it is 
known that such a rule does not exist [1|, once could ask 
which response curves are possible in CA rules? We propose 
to approach this problem from an opposite direction: given 
the CA rule, what can we say about its response curve? It 
turns out that in surprisingly many cases, the response curve 
can be calculated exactly, providing that preimage sets of 
finite strings under the CA rule exhibit recognizable patters. 



2. Definitions 

In what follows, we will be concerned with what we 
call two-dimensional elementary cellular automata, which 
have a local function depending on the central site, its 
right neighbour, and its top neighbour, and which allow two 
states only, and 1. We will say that these are rules with 
"L-neighbourhood", since the neighbourhood has the shape 
of the letter L. Such three-input binary local rules can be 
considered the simplest "truly" two-dimensional CA rules, 
hence the name "elementary". 

Before we define such rules formally, we will first intro- 
duce the concept of triangular blocks, defined as regions of 
2D configurations in the shape of isosceles right triangles. 
The set of triangular blocks of size r, denoted %, is the set 
consisting of elements 

h,r 

. . (1) 

h,i ■ ■ ■ b r ,u 

where each 6j j G Q. The set of eight blocks in T 2 will be 
called basic blocks. 

We may define the local mapping, or local rule, of a 2D 
CA with L-neighbourhood as g : T2 —> G- The local mapping 
g has a corresponding global mapping, G : G^' — > G^ such 
that (G(s)) id = g ( s ^+ 1 , 1+1J ) , for any i, j G TL, s G G^ ■ 

The block evolution operator g : T r —> T r -i will be de- 
fined as a function which transforms triangular block (fTh into 
another block, c G where c,-j = g (^ b ^ J+1 b +1 J G G 

for i G {1, . . . , r — 1} and j G {1, . . . , r — i). We denote 
g" : T r +n — > T r to be the operator obtained by composing 
g with itself n-times. 

Occasionally, we will need to define a distance between 
two configurations. One can show easily that for s,t G £7 Z 
and i,j G Z, the following satisfies all axioms of a metric: 



d(s,t)=l 1 + ( max {KUil} : s i,j ^*t,j)) 

[0 if s = t 

For 2D CA with L-neighbourhood, we adapt the numbering 
system used in J2|. A local rule g is assigned a Wolfram 
number W as follows 

W(g)= J2 9(Za 2 )2' la «+ 2 ^ +a \ (2) 

ao,ai,a 2 G{0,l} 



We note that, as in the case of radius- 1 ID CA, there are 
256 possible elementary 2D CA. Many of these rules are 
related to each other by the group of 4 transformations 
D\ x £?2, where D\ is the dihedral group with a single 
refiectional symmetry and 52 denotes all permutations of 
the elements in {0, 1}. Among each class of four (not 
necessarily distinct) rules, we choose one representative with 
the smallest Wolfram number. We denote this rule to be 
a minimal rule. A list of all 88 minimal rules and their 
equivalences can be found in J3|. 

3. Densities of Blocks 

We now consider the density response problem. Suppose 
that we start with an initial configurations in which a certain 
proportion of sites is in state 1. The simplest way to achieve 
this is to set each site to be in state 1 with probability p, 
and with probability 1 — p, doing it independently for all 
sites. This means that the probability of randomly selected 
site to be in state 1 is p. Suppose that we apply n iterates of 
some CA rule to such configuration. What is the probability 
that in the resulting configuration, the state of a randomly 
selected site is 1? 

In order to formulate this problem more precisely, we will 
use the concept of probability measure, similarly as done in 
[4 1, for one-dimensional CA. 

Given a block b 6 7^-, we define a cylinder set given by b, 
Cij(b), as the set of all configurations in which block b is 
fixed and placed at coordinate aligned at the lower- left 
element of b. We define a measure of such as cylinder set, 
fj, [Cij (b)], to be the probability of occurrence of block b 
placed as above. If the measure is translationally invariant 
we may drop the indices For p g [0, 1], the Bernoulli 
measure is a measure where all sites are independently set 
to 1 with probability p, and to with probability 1 — p. In 
such case, we have 

p p [C(b)]=pi(l-p)( r2+r V 2 -J, (3) 

where j is a number of cells in state 1 in b. 

We now consider the action of the global mapping G on 
the measure of a cylinder set given by block b, which yields 

(G Mp ) [(7(6)] = p p [G- 1 (C(b))] . (4) 

Considering instead n iterations of G, we obtain 

(G> p ) [C(b)] = p p [G~ n (C(b))] . (5) 

If we let g~™(6) be the set of all n-step preimages of block 
b, that is, the set of all blocks a such that g"(a) = b, then 
we can write 

p p [G- n (C{b))} = »p M ■ <® 

a£g-"(h) 

Using the notation P n (b) = (G n p p ) [C(b)}, we write ^ as 
P„(6)= Y P o( fl )- W 

QGg-"(fc) 



If b = 1, and if the initial measure is Bernoulli, then in 
the above formula each Pq(cl) depends only on p, where 
p = Pq(1). -Pn(l) can then be interpreted as the density 
of Is in the configuration obtained by iterating the CA rule 
n times starting from disordered initial configurations with 
density of ones equal to p. 

Plot of P n (l) versus p will be called a response curve for 
each elementary 2D CA. In the special case when p = 1/2, 
the probability of any block of a given size is equally likely 
and (|7]i can be expressed as 

P n (b) = 2- < - r+n+1 ^ r+n ^ 2 card [g _n (6)] , (8) 

where card [g~™(6)] denotes the number of elements in the 
set g~ n (b). If we want to indicate that we consider the 
special case of p = 1/2, we will use the notation P„ (6), at 
the sequence of P„ (6) for n — 0,1,2... will be called 
a response sequence. Finally, we denote P(b) to be the 
asymptotic density of block b, which we obtain by taking 
the limit of P n (b) as n — > oo (if the limit exists). 

4. Theoretical Response Curves 

For 26 minimal rules, we were able to determine an 
explicit response curve formula. In some cases, we found 
that the response curve was independent of n. In other cases, 
the response curve was dependent on n, and then a separate 
formula for the asymptotic density could be obtained. We 
present in detail three examples of each types. In each 
example, we describe the structure of the preimage sets but, 
due to space constraints, we omit direct proofs while noting 
that each case can be proved easily by induction. 

4.1 Rules with Constant Density 

In each of the following examples the formula for the 
response curve has no dependence on n. Therefore, the 
formula for the asymptotic density is the same as the 
response curve. We provide detailed analysis for Rules 0, 
3 and 42 and the remaining results are presented in Table [T] 

Proposition 1: The response curve for Rule is P n (l) = 

0. 

Proof: There are no triangular blocks of any size 
that can be mapped under g l to single block 1. Therefore, 
card [gQ~™(6)] = and we apply (|7ji to obtain our result. ■ 
Proposition 2: The response curve for Rule 42 is 

P n {l)=p{l-p)(l + p). 
Proof: It can be shown by induction that the only 
blocks that map to a single 1 under g" 2 are either blocks 
in T n where b n l — 1, £>n— 1,2 = and all other elements are 
arbitrary, or blocks in T n where b n< \ = 6 n -i,2 = 1, &n-i,i = 
and all other elements are arbitrary. Using ([3]l and Q, we 
conclude that P n (l) = /°)+/? 2 (l~ p)> which simplifies 
to the desired result. An experimental curve confirming this 
result is presented in Figure [id] ■ 



(a) Rule 2 



(b) Rule 3 (c) Rule 10 (d) Rule 42 

Fig. 1: Experimental Response Curves 



(e) Rule 138 



Table 1: Constant Density Rules 



Rules 




P^ S) (1) 











2 


p(i - pY 


1/8 


3 


(i-p) a 
p(i + p-p 2 ) 


5/8 (n odd) 
1/4 (n even) 


4 


p(i - rt" 


1/8 


5 


(i-p) a 
p(i + p-p 2 ) 


5/8 (n odd) 
1/4 (n even) 


10 


p(i-p) 


1/4 


12 


p(i-p) 


1/4 


51 


1 — p (rt odd) 
p (n even) 


1/2 


34 


P(l " P) 


1/4 


42 


p(l-p)(l + p) 


3/8 


51 


1 — p (n odd) 
p (n even) 


1/2 


76 


p(l-p)(l + p) 


3/8 


170 


P 


1/2 


204 


P 


1/2 



Proposition 3: The response curve for Rule 3 is 



^n(l) = 



(1 — p) 2 if n even, 

p(l + p — p 2 ) if n odd. 
Proof: Since Rule 3 has period-2 behaviour, we must 
consider cases when n is odd and when n is even. When n is 
odd, the only blocks that map to a single 1 under n-iterations 
of g 3 are blocks in % where &i,( n+ i)/ 2 = &i,(«+3)/2 = 
and all other elements are arbitrary. 

When n is odd, the only blocks that map to 1 under g 3 l 
are blocks in T n where &i ( n +2)/2 = and all other elements 
are arbitrary, or blocks in T n where b l n / 2 — &i,( n +4)/2 ~ 
l>^i,(n+2)/2 = and all other elements are arbitrary. 

Using Q we can determine the initial probability of 
occurrence of blocks of either type and using (j7]l obtain a 
formula which simplifies to the desired result. Note that in 
the special case when p — 1/2, we conclude that P„ (1) 
equals 1/4 when n is even, and equals 5/8 when n is odd. 
We averaged our experimental results over an even number 
of time steps, thus Figure [lb] has the form 

Pn(l) = l [(1-P) 2 +P(l + P-P 2 )}. 



4.2 Rules with Decaying Density 

In each of the following examples the formula for the 
response curve is dependent on n, and thus we can also de- 
termine an asymptotic density formula. We provide detailed 
analysis for Rules 32, 128 and 138, while the remaining 
results are presented in Table [2] 

Proposition 4: The response curve for Rule 128 is 

P„(l)=p(" 2 + 3 " +2 )/ 2 . 
Proof: The only block mapping to a single 1 under 
g™28 is me block consisting entirely of ones. We use ([3]) to 
find the initial probability of this block and |7]) produces our 
result. Thus, the asymptotic density under Rule 128 is 



P(l) - lim P„(l) 



if p * 1, 
if p = 1. 



Proposition 5: The response curve for Rule 32 is 

P n {l)=p n+1 {l-p) n . 
Proof: Under rule 32, the only blocks that map to a 
single 1 under gg 2 are of the form 

1 

o ■•. 



* ■ ■ ■ * 1 ■ 

Using ([3]), we can determine the initial probability of 
occurrence of blocks of this type and using |7]), we obtain 
our result. The asymptotic density is thus P(l) = 0. 

■ 

Proposition 6: The response curve for Rule 138 is 

P n l ) = P J/ . 

+ 1 

Proof: The only blocks that map to a single 1 under 
n-iterations of g 138 are comprised entirely of arbitrary 
elements in the top n — 1 rows, and have their lower two 
rows, where i ranges from 1 to n + 1, of the form 

★ ■••★01 ■■•11 
■■■Ill 



Table 2: Density Decaying Rules 



Rules 


Pn(l) 


P(l) 


8 







32 


p n+L (i - p) n 





40 


2'V n+1 (l - p) n 





72 


2 n p n+1 (l - p) n 





128 


p (n 2 +3n+2)/2 


if p ^ 1 

1 if p = 1 


130 


see [5] for complete analysis 


132 


can be derived from Rule 130 


136 


p n + l 


if p # 1 

1 if p = 1 


138 


p 2 " + 2 +P 
P+l 


1 if p = 1 


140 


P 2 " + 2 +P 
P+l 


1 if p = 1 


160 


p n+l 


if p ^ 1 

1 if p = 1 


162 


P 2 "+ 2 +P 
P+l 


1 if p = 1 



Using Q we can determine the initial probability of occur- 
rence of blocks for each possible value of i. Summing over 
all i and using (FTl, we conclude that 



Pn(l) = P 



2n+l 



2n+l 




P 



which simplifies to our desired result. 

Again, we can find the asymptotic density as 



P(l) = lim P n (l) = { 1+ p 

n— >oo \ 



if P t 1, 

if p = 1. 



This result is confirmed by the experimental curve in Figure 
[Te| Note that while the response curve is continuous, the 
asymptotic density has a discontinuity at p = 1, correspond- 
ing to an initial condition consisting entirely of ones. 

5. Theoretical Response Sequences 

In some cases we were unable to determine an explicit 
expression for the response curve of a given rule, but we 
were able to derive an explicit formula for card [g~ n (l)], and 
thus use ([8]) to obtain a response sequence. For 21 additional 
rules, we were able to either prove or conjecture a response 
sequence. We first consider the class of surjective rules. 

5.1 Surjective Rules 

Sites belonging to the L-shaped neighbourhood ( al[l ai,o ) 
will be identified by their indices as (0, 1), (0, 0), and (1, 0). 
Similarly as in J6|, a local function g will be called permutive 
with respect to the (0, 1) site if for any choice of y, z G Q 
the function x — > g(y z ) is one-to-one. Permutivity with 
respect to the central site (0, 0) or the right neighbour (0, 1) 
is defined similarly. We now find a response sequence for 
rules permutive with respect to site (0, 0). 



Proposition 7: The response sequence for Rules 15, 30, 
45, 51, 54, 57, 60, 90, 105, 106, 108, 150, 154, 156, 170 
and 204 is P^\l) = 1/2. 

Proof: There are 16 rules permutive with respect to the 
centre site, of which the following 9 are minimal: 51, 54, 57, 
60, 105, 108, 150, 156 and 204. If a rule is permutive with 
respect to (0, 0), then there must exist numbers xq, . . . , X3 G 
{0, 1} such that the local function takes the form 



l{ a° a 2 ) — < 



if ( a? ct 2 ) € { x i x ± 1 5 x 2 i 13 1 } 

5 

1 if foo ~i fz r 1 1 1 

1 11 V ai a 2 ) C \ ^ Q ; xi 1 > x 2 ! x 3 1 J 



(9) 

where xi denotes 1 — xi. Assuming the above form of g, let 
us consider an arbitrary block b G T n - We will now show 
how to construct all preimages of b under g. First of all, we 
claim that blocks c e T" +1 of the form 



Oil 



c i, 



Cl,l 



(10) 



are the only preimages of b, where each a,; (1 < i < n + 1) 
is an arbitrary value in {0,1}, and values of Cjj € {0,1} 
can be determined by an iterative algorithm. 

To see that this is indeed true, we now present an algo- 
rithm with which we can construct all possible preimages: 

1) Starting from bi n , we wish to find all neighbourhoods 
{al a 2 } such that g( l° ia2 ) = b ltn . The structure 
of the local mapping gives us four possible such 
neighbourhoods {ala 2 } = { c?°„ ai }, where 

Cl,n = (1 — i , l,n)l2ai+a 2 + &l,n(l _ X2a 1 +a 2 ), 

and the values of a\ and ct2 are arbitrarily selected. 
We now repeat step 2 for all values of i € {2, . . . , n}. 

2) Since 6; in -;+i is given and ot. L has been freely chosen 
in the previous iteration, we wish to know all neigh- 
bourhoods { H a2 }, such that g( "° Q2 ) = 6 4 , n _ <+1 . 
The structure of the local mapping gives two possible 
neighbourhoods {a% 2 } = { Cl ,„ Q * i+1 a i+1 }, where 

Ci,n-i+l = (1 — bi. n -i+i)x2cti+a i+1 + 
+ &i,n-i+l(l ~ s 2Qi+Q i+ i)j 

and ai+i is another arbitrarily selected value. 
We now construct the rest of the preimage and show 
that all other values are uniquely determined based 
on each choice of the a values in the top diagonal. 
For all values of j G {1, . . . , n — i} and then for all 
i G {1, . . . , n — j}, we repeat step 3 as follows. 

3) Since j+i is fixed, we wish to know all 
neighbourhoods { l\ a2 }, such that g( a% 2 ) = 



Since c, „_j_,- +2 and a 



A-\-l.n- 



-3+1 



fixed in a previous iteration, the structure of the local 



(a) Rule 2 



(b) Rule 108 (c) Rule 154 (d) Rule 168 

Fig. 2: Experimental Response Curves 



(e) Rule 172 



mapping tells us that our neighbourhood must have 
the form { % a , 2 } = { ^Zt^l c 1+1 ,„_,_ J+1 }, where 

and i' = 2c iin _ i _ i+2 + Cj + i,„_i_j + i. Note that no 
new arbitrary parameter appears here, thus the neigh- 
bourhood is determined uniquely. 

The only arbitrary values in the preimage are the (n + 
1) values of on on the main diagonal. Therefore, we know 
that there are exactly 2 ,i+1 preimages for a given b G T™. 
Therefore, we can see that card [g -7l (l)j = 2 (n2+3n+2 '>/ 2 . 
Now, using (jij, we conclude that P„ (1) = 1/2 for all n. 
Considering rules permutive with respect to the other two 
sites, we conclude that also Rules 15, 30, 45, 90, 106, 154 
and 170 possess a response sequence Pn(l) = 1/2. ■ 

It turns out that these rules are the class of surjective 2D 
CA with L-neighbourhoods. In one dimension, it is known 
that rules permutive with respect to one of the variables 
located at the left or the right end of the neighbourhood 
are surjective, as proved in [6|. Recently, this result has 
been generalized to two dimensions by Dennunzio and 
Formenti 0, who demonstrated that any rule with Moore 
neighbourhood (of any radius) which is permutive with 
respect to one of the corner sites is surjective. We now show 
how one can prove a similar result specifically for the L- 
shaped neighbourhood, adapting the idea in [8| to 2D CA. 

Proposition 8: If the local mapping of an elementary 2D 
CA with L-neighbourhood is permutive with respect to any 
site, then the corresponding global mapping is surjective. 

Proof: From Proposition [7] we know that for any 
permutive rule and all b G T n ,(n > 1), card[g _1 (i>)] = 
2 n+1 . Consider any infinite configuration, t G G z ■ Define 
for all n > 1, the set, S n = {s e G z : g(s [n+1 ]) = t [n] }, 
where s\ n+ x\ denotes a block of size n contained in an 
infinite configuration s E G z and placed at (0,0). Our 
assumption guarantees that all S n are non-empty for n > 1. 
We also know that S n +i Q S n . We consider the complement 
of S n , the set S n = {s £ G z : g(s[„+i]) ^ t [n ]}, to show 
that S n is a clopen set. 



We first show that S n is open. Let s € S n C {0, 1} Z be 
an arbitrary configuration. For all e > 0, we choose k E Z, 
where k > n, such that < e. We now pick an infinite 
configuration s' 6 {0, 1} Z such that d(s, s') = k , 1 +1 , where 
k* > k. Since s £ S n , we know that s' € S n , and 



Thus, S n is open. Similar analysis shows that S n must also 
be open, and thus S n is a clopen set. By the Nested Set 
Theorem [9|, there must exist s s Q z , such that F(s) = t. 

m 

To conclude that these are the only surjective rules, we 
use the reverse direction of the Balance Theorem. 

Proposition 9: If a elementary 2D CA with L- 
neighbourhood is surjective, then for all n > 1 and 
all blocks b G T\ cardfe" 1 ^)] = 2" +1 . 

Proof: The Balance Theorem was proved in ID in O 
and in 2D in ifTUl . A version of the proof specifically tailored 
for the L-neighbourhood is to be reported elsewhere 0. ■ 

For all other elementary rules, we performed a computer- 
ized search and found blocks for which no preimages exist. 
By Proposition [9| these rules must be non-surjective. 

5.2 Conjectured Response Sequences 

To find response sequences for the remaining rules, we 
performed an exhaustive search through all potential preim- 
ages for each rule. For the L-neighbourhood, the number of 
potential preimages is 2(™"+ 3n+2 )/ 2 , which makes searches 
for large n impossible. We performed our searches using 
the Shared Hierarchical Academic Research Computing Net- 
work (SHARCNET) and we were able to obtain cardinalities 
of preimage sets to level n = 7. We then attempted to 
conjecture a formula for the sequence using the first six 
terms, and checked the conjecture with the seventh term. 

Rules 23, 27, 29, 43, 46, 58, 77, 78, 142, 172, 178, 184 
each shared the first seven terms of the preimage sequence 
with the surjective rules above, so that for these rules we 
conjecture that P n s \l) = 1/2. For all remaining rules, a 
list of the first seven preimage cardinalities is available upon 
request. 



6. Experimental Response Curves 

For those rules for which an explicit response curve 
formula could not be derived, we were able to perform com- 
puter simulations to obtain experimental response curves. 
We start with a square configuration of 250000 elements 
and we iterate 1000/p(l — p) times when p g (0,1) and 
100000 times otherwise, with periodic boundary conditions 
and averaging density over the last 10 time steps and over 10 
iterations from different initial conditions. Examples of such 
experimental curves are presented in Figures [T] and [2] We 
note in passing that one of the examples shown in Figure |2j 
namely Rule 168, exhibits response curve resembling "phase 
transition", that is, discontinuity of the derivative. None of 
ID elementary rules exhibits such behaviour. 

7. Rule emulation 

We will now briefly turn our attention to dynamics of 2D 
rules. When one prints sample spatiotemporal diagrams of 
2D rules with L-shaped neighbourhood (not shown here for 
the lack of space), one can easily observe that all rules for 
which density response curves can be calculated theoretically 
exhibit somewhat "simple" dynamics. A convenient was to 
describe this "simplicity" is to say that after a few iterations 
these rules essentially behave like identity or shift. In order 
to formalize this statement, we need to introduce the concept 
of emulation, first finite and then asymptotic. 

7.1 Finite Rule Emulation 

We say that Rule X emulates Rule Y at level n if, 

Ex + \b)=g Y (g x (b)). (11) 

for any block 6 £ B n+2 . We will demonstrate this with an 
example. Consider Rule 76, with a local rule given by 

976 (1 z) = {l-x)y(l-z)+(l-x)yz+xy(l-z) = y(l-xz). 

We now compose gjg with itself as follows 

&Q>)=gn (gw(g g B5 )) 

= x 3 (l - xix 4 ) (1 - £1(1 - x X 2 )x4(l - x 2 x 5 )) 

= £3(1 - X1X4) = .g 20 4 (g76(&)) ) 

where we have used the fact that when x G {0, 1}, we 
know that x 2 = x. We therefore conclude that Rule 76 
emulates identity at level 1. We checked all 88 x 87 pairs 
of distinct elementary rules for finite rule emulation. In 
Figure [3] we show all level 1 emulation relations between 
all minimal elementary 2D rules with L-neighbourhood as 
directed graphs in which an arrow travels from X to Y if 
and only if Rule X emulates Rule Y at level 1. In Figure 
[3a] are all rules which finitely emulate the identity Rule 204. 



In Figure 3b are all rules which finitely emulate the left 



shift Rule 170. Finally, in Figure 3c are another class of 
interrelated emulation rules. In addition to the rules in the 
graph, we also discovered that rules 6, 14, 18 and 50 emulate 
rules 134, 142, 146 and 178 respectively. 



7.2 Asymptotic Rule Emulation 

In ifTTI . the author defined the following metric to describe 
the distance between two elementary ID cellular automata 
rules. We adapt this and define the following metric to 
describe the distance between two elementary 2D cellular 
automata rules with L-neighbourhood 

d(/„g)=2(- fc2 - 3fc - 2 )/ 2 ^ . (12) 

beB k 

We say that Rule f asymptotically emulates Rule g if 

lim d(f n+ \f n og) =0. (13) 

n— > 00 

We now derive a useful equation with which we can calculate 
the distance between two rules at a given level-n. First, we 
define the following function for any block b E B, 

(f®g)(b) = f(b) + g(b) mod 2, 

which outputs 1 if and only if f(b) 7^ 5(6). Thus, we can 
use this function to count the number of blocks on which 
local mappings / and g differ. Adapting Proposition 3 from 
IfTTI . we obtain the following proposition (proof in fTD ). 

Proposition 10: If /, g are 2D local L-neighbourhood 
mappings, A = (f g)~ x (1), and A n = f~ n (A ), then 

card [A n ] 



d(f n+1 J n og) = 



(14) 



2(« 2 +5n+6)/2 - 

We demonstrate this procedure with an example. 

Proposition 11: 2D CA Rule 160 asymptotically emu- 
lates the identity rule. 

Proof: If we consider the local mappings for both 
Rules 160 and 204 we see that the set of blocks on 
which the rules output differ is^4o = {io,ii,oi'io} 



To use Proposition 10 we must find the set A n in gen- 
eral and thus we must know the sequence of preimages 
for these particular four basic blocks. We found the first 
five terms of these sequences and conjectured patterns are 
B (2"+ 2 -3) ,B,B (2" +1 - 1) and B, respectively, where 
B = 2(" +™)/ 2 . From equation (14i, we determine that 



d (§160 . §160 §204) =3-2 - 4 

Therefore, since the limit of this expression goes to 0, we 
conclude that Rule 160 emulates identity asymptotically. ■ 
Table [3] shows all known results of rules emulating shift 
or identity. We can now state our observation expressed at 
the beginning of this section using the concept of emulation: 
all rules included in Table s^and^emulate identity or shift 
either in a final number of steps or asymptotically. 

8. Further Results: Basic Blocks 

We also note that if p = 1/2, it is often possible to 
compute the number of preimages of other blocks. For 
example, for 40 of the 88 minimal rules, we were able to 
find preimage sequences for all eight basic blocks, that is, 
blocks in 72- In each case, it is only necessary to determine 






(a) Rules Emulating Identity (b) Rules Emulating Shift 

Fig. 3: Finite Emulation Relations 



(c) Other Finite Emulations 



Table 3: Asymptotic Emulation 



Rule / 


d(/ n+ \/204 °/ n ) 


pW(l) 


8 


3 . 2-' Zn -~ 6 





32 


5 ■ 2-' 2n -' i 





40 







72 


4— n— a 





128 


2(-n^-3n-2)/2 _ 2(-i 2 -5n-6)/2 





132 


2 (-n i! -3n-4)/2 


~ 0.179 


136 


2 -n-2 





140 


2~ 2n — 3 


1/3 


160 


3 • 2 _n ~ 2 - 4-"- 1 





Rule / 


d(f n+1 , f 170 of n ) 


PW(1) 


130 


2 (-n 2 -3n-4)/2 


~ 0.179 


138 


2~ ^ro — 3 


1/3 


162 


2~ 2n— 3 


1/3 



preimage sequences for 5 of the 8 blocks, then we may use 
Kolmogorov consistency conditions |[T2l to determine the 
remaining three. In some cases, these formulas are rather 
striking, such as in the case of rule 1 30, reported in detail in 
0, or rule 172, for which we make the following conjecture. 

Conjecture 1: Under 2D CA Rule 172 the preimage se- 
quences of basic blocks are given by 



card [g 



2(n 2 +5n)/2 ^ Ck 



k=0 



2 (n*+5n)/2 ( 2 - £ & 



k=0 



if be B x , 
if b € B 2 , 



where Ck denotes the k-th Catalan number and 

D f 1 IT. TD f 1 11 

D \ — t ' 11)00' 11/' 2 I 1 J 10:01' 10J- 

Work on a proof of this result is ongoing and will be reported 
elsewhere. 

9. Conclusions and future work 

We demonstrated that response curves are calculable for 
simple rules that emulate shift or identity. Response curves 
clearly deserve further study and it is worthwhile to sys- 
tematically study them for other CA rules. However, due to 



rapidly increasing preimage size, this won't be an easy task 
for larger neighbourhoods. One would need a more efficient 
way to construct the set of preimages of a given block, 
as simple brute force search becomes computationally too 
expensive. We also hope that rigorous results can be obtained 
for rules with somewhat more complicated dynamics. 
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